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Opening a sizable band gap without degrading its high carrier mobiUty is as vital for silicene as for graphene 
to its application as a high-performance field effect transistor (FET). Our density functional theory 
calculations predict that a band gap is opened in silicene by single-side adsorption of alkali atom as a result of 
sublattice or bond symmetry breaking. The band gap size is controllable by changing the adsorption 
coverage, with an impressive maximum band gap up to 0.50 eV. The ab initio quantum transport 
simulation of a bottom-gated FET based on a sodium-covered silicene reveals a transport gap, which is 
consistent with the band gap, and the resulting on/off current ratio is up to 10^. Therefore, a way is paved for 
silicene as the channel of a high-performance FET. 
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Owing to its extremely high carrier mobility, graphene has attracted enormous attention since its discovery 
in 2004 \ However, pure graphene is not suitable for effective field effect transistor (FET) operating at 
room temperature due to its zero band gap. Opening a sizable band gap without degrading the electronic 
properties is critical for the application of graphene in nanoelectronics and is probably one of the most important 
and tantalizing research topics of the graphene community. Any successor to silicon metal-oxide-semiconductor 
FET (MOSFET) that is used in complementary MOS-like logic must have a sizeable band gap of 0.4 eV or more 
and an on-off current ratio between 10^ and 10^ ^. Unfortunately, the opened band gap without degrading the 
electronic properties in graphene (via a vertical electric field, molecular surface adsorption, or hexagonal boron 
nitride (h-BN) as a substrate) is usually smaller than 0.34 eV^"^, and the observed on-off current ratio at room 
temperature is no more than 100^. Silicene, graphene analogue for silicon, is predicted to be a zero band gap 
semiconductor with a Dirac cone as observed in graphene^°'^\ The synthesis of silicene on the Ag surface^^"^^ or 
zirconium diboride thin film^^ has already been reported in several experiments. Recently, the linear band 
dispersion near the Fermi level has been observed in a silicene sample on Ag(lll) surface by angle- resolved 
photoemission spectroscopy (ARPES), unambiguously verifying the existence of silicene^^. With its massless 
Dirac Fermions character, extremely high carrier mobility is expected in silicene as in graphene. Therefore, 
silicene is a promising material for high-speed switching devices. Apparently, pure silicene cannot directly be 
used for effective FETs either because of the same cause as graphene^^. 

Although covalent functionalization^°"^^ or epitaxial growth in the form of nanoribbon^^'^^ can open a sizable 
band gap in silicene, the electronic properties are destructed, and the carrier mobility is inevitably decreased to a 
large extent, as happened in graphene^^'^'*. Therefore, opening a sizable and tunable band gap of silicene without 
degrading the electronic properties is equally highly desired for high-performance silicene FET devices. It is 
proposed that a vertical electric field can open a tunable band gap in silicene without degrading the electronic 
properties'^. However, the maximum opened band gap by the experimentally accessible electric field is just 
0.13 eV. Moreover, a single-gated device is desirable from the device point of view^^. Very recently, it is reported 
that silicene that epitaxial grown on Ag surface or zirconium diboride thin film can open a band gap of about 0.6 
or 0.3 eV, respectively'^''^. However, Ag is an metal substrate, which is undesired in a FET device, and the 
calculated dispersion relation near the band gap for silicene on zirconium diboride thin film is very flat, which 
will lead to a rather large carrier effective mass and greatly decrease the carrier mobility. Besides, the size of the 
band gap of silicene opened by the interaction with substrate is difficult to tune. 

In this Article, based on the density functional theory (DFT) calculations we propose that the band gap of 
silicene can be opened without degrading the electronic properties by surface adsorption of alkali metal (AM). 
The band gap is tunable by controlling the coverage of AM. Most intriguingly, the maximum band gap is up to 
0.50 eV by adsorption of AM atoms, which meets the band gap requirement of high-performance FETs. 
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Subsequently, a bottom-gated FET device based on a sodium-cov- 
ered silicene monolayer is simulated by the ab initio transport cal- 
culations. A sizable transport gap and a high on/off current ratio 
(~ 10^) are available in this device, suggestive of the great potential 
of silicene as the channel of a high-performance FET. 

Results 

The supercell model consists of an m X m primitive silicene cell (m = 
1, V3, 2, 3, and 4) and an AM atom (AM = Li, Na, K, Rb, and Cs) on 
the top. The AM- covered silicene is denoted as AMSin, where n is the 
number of the silicon atom per supercell. We define the coverage AT of 
the AM on silicene as the ratio of the number between the AM and Si 
atoms per supercell, and one has N = I In. The in-plane lattice con- 
stant of silicene is taken from that of previous works, a = b = 3.866 
^10,26,27 AiQj^g ^j^g z-direction perpendicular to the silicene plane, the 
lattice constant of 30 A is used. 

Geometry and stability. Of all the examined adsorption sites, the 
AM atoms favor sitting above the hexagonal ring center of silicene 
and have the highest sixfold coordination (Fig. 1). Such a sixfold 
coordination configuration of the AM atom has been reported to 
be favorable in a AM-covered graphene^^. A metastable state, in 
which the AM atoms all sit above the bottom Si atom, is available 
(Supplementary Fig. SI). The energy differences between these two 
states increase from 0.06 to 0.21 eV with the decreasing coverage, 
which are larger than the room temperature (~ 0.026 eV) and thus 
ensure the thermodynamically stability and reliable operation if the 
AM-covered silicene is fabricated as the FET channel. 

As shown in Supplementary Fig. S2, the silicene buckling (do) 
increases from 0.46 to 0.49 ~ 0.77 A after the AM adsorption, weakly 
dependent on the AM species and the coverage. On the contary, the 
distance of the AM atom to the silicene top-surface (di) singnifi- 
cantly ranges from 1.3 to 2.8 A with the increasing atomic radius 
of AM, but nearly independent of the coverage N. 

We define the adsorption energy of the AM atom on silicene as 

Esi=Esi-\-EAM—Esi + AM, (1) 

where Esi, £am and £si+AM are the relaxed energy for silicene, the 
isolated AM atom, and the combined system, respectively. As shown 
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Figure 1 | Structures of the Na-covered m X m silicene supercell. (a-e) 
The adsorption coverage is N = 3.1%, 5.6%, 12.5%, 16.7%, and 50.0%, 
respectively, for m = 4, 3, 2, ^/3, and 1. The rhombi plotted in black 
line represent the supercells at different coverages, (f) Top and side views of 
an Na-covered silicene supercell over the h-BN sheet with N = 12.5%. do: 
the silicene buckling; dii the distance between Na and the top -surface 
silicene; <i2- the distance between h-BN and the bottom-surface of 
silicene. 



in Fig. 2a, the E^ values are positive with magnitudes of 0.57 ~ 
2.92 eV per AM atom, showing that the AM one-side adsorption 
on silicene layer is always strongly or moderately exothermic. E^ 
tends to reduce with the increasing N. One cause is that the electro- 
static attraction between the AM atom and silicene is reduced with 
the increasing N, as a result of the reduced charge transfer Q from the 
AM cation and silicene (Fig. 2b). Another cause is the repulsion 
between the AM cations and silicene increases with the increasing 
N. We notice that although the transferred Q from Li to silicene is 
apparently smaller than those from other AM atoms, £a of Li on 
silicene is comparable with those of other AM atoms on silicene. 
This suggests an evident 'covalent component' of Li-Si interaction 
besides the ionic contribution^^. 

Electronic structure. Given the same coverage, silicene covered by 
different AM species share similar electronic structures. We select 
Na-covered silicene as a representative and show its electronic 
structure at different coverages in Fig. 3. Even at the smallest 
coverage of N = 3.1%, a small direct band gap of 0.04 eV is 
opened between the n and 7i* states while the Dirac cone shape is 
almost preserved (Fig. 3a). Due to the heavy electron doping, the 
Fermi level (Ef) is elevated by 0.50 eV with respect to the valence 
band maximum (VBM). The band gap at the Dirac point is increased 
monotonically to 0.08, 0.17, 0.30, and 0.45 eV atiV = 5.6, 12.5, 16.7, 
and 50%, respectively (Fig. 3b-3e). Notably, both of the two original 
Dirac points and K'^ in silicene are folded to the T point of the 
reduced first Brillouin zone (BZ) at iV = 5.6% {^/3x^/3 supercell) 
and 16.7% (3X3 supercell), as shown in Fig. 3b and 3d, respectively. 
The opened maximum band gap has already met the requirement of 
0.4 eV for practical FETs operation. AtN = 12.5 and 50%, there is a 
band across £y around the T point, which mainly originates from the 
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Figure 2 | (a) Adsorption energy, (b) Mulliken charge transferred from 
the adatom to silicene, (c) band gap, and (d) effective mass of electrons of 
the AM-covered silicene as a function of the coverage N. Lines connecting 
symbols are guides to the eye. 
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Figure 3 | Band structures of the Na-covered silicene at different 
coverages. The Dirac point in (b) and (d) is folded to the F point due to the 
reduction of the first Brillouin zone. The red dashed-line in (c) stands for 
the band structure of the Na-covered silicene with an h-BN buffer layer. 
Inset of (e): Orbital of the conduction band bottom at the T point. The 
Fermi level is set to zero. 

broaden outmost 3s valence orbital of the adsorbed Na atoms (The 
inset of Fig. 3e). 

Fig. 2c shows the band gap change with the coverage N of the five 
kinds of the AM- covered silicene. Given a coverage the band gaps 
induced by different AM species adsorption are similar except for the 
LiSig and LiSie monolayers, the latter of which has an extraordinary 
large band gap of 0.46 eV, compared with a value of about 0.2 eV for 
other AM-covered silicene at this coverage. The band gaps of all the 
five kinds of AM-covered silicene monolayers increase generally with 
the increasing coverage, with values ranging from 0.04 to 0.50 eV. 

It has been established that a vertical electric filed can induce a 
band gap in silicene because it breaks the sublattice symmetry in 
sihcene^^. The AM adsorption leads to a charge transfer between 
AM and silicene, which builds a perpendicular electric field in sili- 
cene. This perpendicular electric field breaks the sublattice symmetry 
and thus opens a band gap in silicene. The asymmetry between the 



two sublattices increases with the increasing coverage, and the band 
gap is correspondingly increased. We should notice a similar phe- 
nomenon and band gap opening mechanism in bilayer graphene. 
Both experiments^'^'^^ and theoretical calculations^'^^ reveal that a 
band gap can be opened in bilayer graphene via single-side adsorp- 
tion of metal atoms or acceptor/donor molecules because a vertical 
electric field is built in bilayer graphene due to a charge transfer 
occurs between the adsorbate and bilayer graphene and it breaks 
the inversion symmetry of bilayer graphene. The size of the band 
gap can be controlled by changing the adsorption concentration. 
Besides the breaking of sublattice symmetry, we point out that the 
breaking of the bond symmetry also contributes to the band gap 
opening in the AM-covered silicene. Let us recall its single layer 
graphene counterpart. When adsorbed on single layer graphene, 
the Li atoms are also located above the hexagon center of single layer 
graphene (notated as LiCn), just like the case happened in silicene. If 
and only if in some certain coverages (n = 3x, x is an arbitrary 
integer) and the corresponding supercell take forms likeV^ x \/3, 
3X3, etc. (The lattice constant is their integral multiple), the bond 
symmetry is broken, resulting in a Kekule ordering. Finally, a band 
gap is created in the LiC^^ monolayer though the sublattice symmetry 
of single layer graphene is well preserved^\ 

In an AMSin monolayer, when n = 3x, the second band gap 
opening mechanism also works. It is convenient to interpret this 
mechanism using the tight-binding model with the Hamiltonian 



(2) 



where c] and q are the fermionic creation and annihilation operators 
on site /, £/ the on-site energy, and t^^ the hopping amplitudes. 
Following this model, we can divide the effects of AM-adsorption 
into two parts: 

1) Breaking of sublattice symmetry. The on-site energies 8^ are 
equal in the pure silicene. Due to the charge transfer between 
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Figure 4 | Kekule structures of AM-covered silicene. (a-b) Electrostatic 
potential distributions in Li and Cs-covered silicene with N = 16.7%. The 
rhombi plotted in black line represent the supercells. (c) Schemetic model 
of the AMSig supercell. ti and t2 are hopping parameters. Top and bottom 
sihcene sublattices are colored yellow and blue, respectively. The numbers 
are labeling of the atom sites, (d) The fc-space of the AMSi6 supercell. The 
Brillouin zones of the pure silicene and the AMSi6 monolayer are colored 
purple and green respectively. The K and FC Dirac points of silicene can 
reach each other via Gi , G2 (the unit vectors of the reciprocal lattice of 
the AMSig monolayer) or their linear combination. 
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AM and silicene, a perpendicular electric field is built in silicene. 
Because there is buckling in silicene, the Coulomb fields felt by 
the top and bottom sublattices are different, resulting to a split of 
£i. We notate the on-site energy of the two sublattices as Si and 82, 
respectively. The eigen function can be approximately expressed 
as a linear combination of the normalized orbitals of Si atoms. 
We consider the most simple case of n = 2 of the AMSin mono- 
layer, and the Hamiltonian is given by 



Sodium 



Silicene 



H = 



- ikaj 



£2 



(3) 



where the vectors (Ti connect the two sublattices of silicene 

(Fig. 4c). In the Dirac points {t 1 = 0)' ^he opened band 
gap can be written as 

A=|£i-£2| (4) 

Only \i n 7^ 3x, the band gap always can be expressed as the 
function of the on-site energy s^. The £/ discrepancy between the 
two sublattices is also responsible for the band gap opening of 
silicene under an external vertical electric field^^. Apparently, the 
AM adsorption on the one side of silicene is equivalent to 
application of a vertical electric field E^{x, y) with periodic dis- 
tribution on the silicene plane. 

Breaking of bond symmetry. Apart from the diagonal terms, 
there are off- diagonal terms in the Hamiltonian that affect the 
electronic structure. When n = 3x, there is a Kekule modulation 
of the hopping amplitudes t^y In pure silicene, t^^ is a uniform 
value (= t).As shown in Fig. 4a and 4b, the potential of bonds in 
the hexagon with the AM atoms is quite different from those 
without them. As a result, the bond symmetry is disturbed, and 
the nearest-neighbor hopping parameter is split into two differ- 
ent values ti and ^2- Take the AMSie monolayer as an example 
(Fig. 4c). Considering their distance respected to the AM atom, 
we define the on-site energies of the top and bottom Si atoms as 
£1 and £2 respectively. Since the Dirac points of silicene are folded 
to the r point of the supercell, the eigenstates and eigenvectors 
can be simpUfied as: 
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where Ui and bi are the coefficients of the normalized orbitals 
of the top and bottom Si atoms respectively in the linear combi- 
nations of the wave function (see Fig. 4c for the labeling of the 
atom sites). The resulting band gap A = \ti — t2\, which is 
independent of £1 and £2. Therefore, the band gap opening at 
the Dirac point for the AMSisx monolayer is completely attrib- 
uted to breaking of the bond symmetry. 

Back to Fig. 3c, the band gaps of the LiSi6 and LiSiig monolayers 
are found apparently larger than those of other AM-covered mono- 
layers at the same coverage. This difference can also be attributed to 
the stronger bond imbalance by Li adsorption. Apart from the ionic 
interaction between the AM atoms and silicene, the AM atoms have a 
dative covalent contribution, which adds extra imbalance to the 
bonds. The covalent contribution is the most significant for Li and 
decreases from Li to Cs^^. As a result, the bond symmetry is broken 
the most, and the magnitude of \ti — t-^ is the maximized by Li 
adsorption. Comparing Fig. 4a with 4b, we can see that the potential 
discrepancy of the LiSie monolayer between the bonds in a hexagon 
with and without the AM atoms is greater than that of the CsSie 
monolayer. What we want to emphasize here is that the Kekule 




Figure 5 | Schematic model of the FET based on the Na-covered silicene. 

The channel is 113.8 A long, and the electrodes are composed of semi- 
infinite silicene. 

distortion happens only if n = 3x. For instance, no band gap has 
been experimentally observed in Na-covered graphene sample 
(NaCg) grown on SiC substrate^^. 

We can also interpret the "three" dependence of the band gap 
creation in silicene in the /c-space. The AM atoms serve as two-dimen- 
sional periodic potential, and certain distribution pattern of the AM 
atoms imposes various periodic boundary conditions into the system. 
In our case of the AMSisx monolayer, the AM atoms expand the lattice 
period and shrink the period in the corresponding Brillouin zone 
(BZ). As shown in Fig. 4d, the reduced BZ of the AMSie monolayer 
is colored green, and the purple one indicates that of silicene. Both the 
two inequivalent Dirac points of silicene and iCg are mapped on the 
r point by the reciprocal unit vector Gi or G2 in the reduced BZ (An 
extra 71/6 rotation with respect to the original BZ of silicene is also 
included for the AMSie monolayer). There is interaction between 
and iCg, and thus the intervalley hopping matrix element is normally 
nonzero, which opens a band gap in silicene. 

Effective masses of electrons (me) of the AM-covered silicene are 
slightly anisotropic, with a discrepancy less than 5%. The anisotropic 
propagation of charge carriers is related to their chiral nature. An 
similar anisotropic behavior of the Dirac fermions under periodic 
potential in graphene has already been reported by Louie et a?^. The 
me averaged over the different directions is provided in Fig. 3d. The 
me values generally increase with the coverage except Li, with values 
ranging from 0.02 to 0.20 mo (mo is the free electron mass). The values 
of me are more dispersed at AT = 5.6 and 16.7% than dXN = 3.1, 12.5, 
and 50.0% because the band structures near the Dirac point at the 
former two coverages are dominated by the AM- sensitive breaking of 
bond symmetry while those at the latter three coverages originate 
from the AM-insensitive breaking of sublattice symmetry. In light 
of the fact that suspended bilayer graphene has an effective mass of 
me = 0.03 mo and carrier mobility of ^ = 2X10^ cmW•s^^ the 
carrier mobility fi in AM-covered silicene is estimated in the range 
of 2X 10^ ~ 10^ cm^/V»s, assuming that its scattering time t is similar 
to that of suspended graphene. Since the carrier mobility of bilayer 
graphene is not significantly degraded by Al and molecule single-side 
adorption^'^^'^^, the scattering time t in silicene is probably not sig- 
nificantly decreased by the AM adsorption. We have also calculated 
the effective mass of holes (mh) and found the average mh is 1% ~ 9% 
less than the average me, and the asymmetry between electrons and 
holes becomes more obvious with the increasing coverage. 

Transport property. To assess the transport performance of the AM- 
covered silicene, we further simulate an FET out of the Na-covered 
silicene with N = 50%. As shown in Fig. 5, the central region of the 
simulated FET is a NaSi2 monolayer with a length of 113.8 A. The 
electrodes are composed of semi-infinite silicene. The dielectric 
region, with a thickness of 7 A, consists of a Si02 substrate and h- 
BN buffer layers to preserve the high carrier mobility of silicene^^. 
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Since the h-BN and Si02 share similar dielectric constant (g), we use a 
homogenous 8 = 3.9 to the dielectric region in the transport 
calculation. A bottom gate is placed below the dielectric region. 

When Vg = 0 V, the transport gap is located below Ef and about 
0.7 eV in size with the SZ basis set (Fig. 6a), a value about 1.5 times 
the band gap of 0.45 eV obtained with the DNP basis set. The SZ 
basis set that used in the device simulation may be inefficient in 
describing some systems and usually overestimates the transport 
gap^^. AppUcation of a larger SZP basis set results to an improved 
transport gap of 0.40 eV (Supplementary Fig. S5). Within the bias 
window, there is a broad transmission peak, which is responsible for 
the on-state current of 236 |iA/|im. As decreases to — 18 V, £y 
shifts downward and becomes closer to the transport gap (See right 
insets of Fig. 6a). Correspondingly, the peak in the transmission 
window becomes narrower and the current decreases to 18.0 [lA/ 
|im. When Vg= —30 V, Ef is located in the opened gap of the Dirac 
cone, and the transmission probability nearly vanishes within the 
transmission window. As a result, an effective off-state with a tiny 
current of 5.45 X 10"^ jiA/jim is achieved. 

Further support for the successful switch capability is provided by 
the transmission eigenstate at Ef and /c-point = (0, 1/3) compared 
between the on- and off-state, as shown in Fig. 6b. The transmission 
eigenvalue in the on -state is 0.70; correspondingly, the incoming 
wave function is slightly scattered and most of the incoming wave 
reaches the other lead. By contrast, the transmission eigenvalue in the 
off- state is merely 0.01; correspondingly, the incoming wave function 
is almost completely scattered and unable to reach the other lead. 

Fig. 6c shows the transfer characteristics of the NaSi2 monolayer 
FET, typical of n type doping. The on/off current ratio in the Na- 
covered silicene FET is as high as 4 X 10^, which meets the require- 
ment for high-speed logic applications. Such an on/off current ratio 
is impressive when compared with graphene, where the observed on/ 
off current ratio without degrading the electronic properties in gra- 
phene is no more than 100 at room temperature and 2000 at 20 K^"^. 
Very recently, the measured on -off current ratio is greatly improved 
to 10000 at room temperature in graphene tunneling FET^^. 
However, the on- state current, which is one important figure of merit 
of an FET and determines the device operation speed, is greatly 
depressed in graphene tunneling FET due to its tunneling feature^^. 
The calculated subthreshold swing from Fig. 6c is S = 140 mV/dec, 
and the terminal transconductance is — 12.6 |iS/|im, ensuring a 
quick current response to the variation in Vg value. 

Discussion 

The h-BN buffer layer is shown to prevent silicene from undesired 
bonding with the Si02 substrate, keeping the high carrier mobility in 
silicene^^. We estimate the effect of the h-BN buffer layer on the 
electronic structure of a freestanding NaSis monolayer. The addition 
of the h-BN sheet seldom disturbs the geometry of the freestanding 
Na- covered silicene since the relatively weak interaction associated 
with a large distance of — 3.6 A (Fig. If). As indicated by the 
redline in Fig. 3 c, the band structure in the vicinity of the Dirac point 
is marginally perturbed by the h-BN buffer layer. Far away from the 
Dirac point, there are some slight shifts of the energy bands near the 
r point, which are of little influence upon the transport performance. 
The h-BN layers induce slightly more asymmetry to the sublattices of 
the Na-covered silicene, involving a tiny band gap increment of 
0.024 eV. On the other hand, it slightly slows down the slope of 
the bands near the Dirac points, leading to a slight increase of me 
from 4.6 X lO"" to 5.5 X 10"" mo. 

We have carried out a calculation including the spin-orbit coup- 
ling (SOC) effects in the NaSis monolayer. The energy differences 
between the most stable and metastable states are nearly unchanged 
(increased by 0.03%). The SOC effects is predicted to cause a semi- 
metal to semiconductor transition in pure silicene with a direct gap of 
— 15 meV in the Dirac point and tends to be larger with heavier 
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Figure 6 | Device performance of the Na-covered (iV= 50%) silicene FET. 

The bias voltage is set to 0.1 V. (a) Transmission spectra with Vg = 0, — 18 
and —30 V, respectively. The vertical dashed-line indicates the bias 
window. The details in the bias window are provided in the left insets. The 
right insets are sketches illustrating the relative positions of Ef and the gap 
manipulated by Vg. (b) Transmission eigenstates (at £ = Ef and /c-point = 
(0, 1/3)) of the on-state ( Vg = 0 V) and off-state (Vg = -30 V). The 
isovalue is 0.2 a.u. (c) Transfer characteristic in linear (left axis) and 
logarithmic scales (right axis). 

buckling^^'^^. As shown in Supplementary Fig. S3, the band structures 
of a NaSis monolayer with and without the inclusion of the SOC 
effects are nearly the same. The band gap increases by only 1.2 meV 
after the inclusion of the SOC effects. 

Many-body effects may significantly enlarge the band gap in semi- 
conductors, especially in the low dimensional systems^^. We also test 
many-body effects on the opened band gap in AM-covered siUcene. 
Supplementary Fig. S4 shows the band structures of the NaSi2 mono- 
layer with and without the inclusion of the many-body effects. The 
71* state around the Dirac points shifts upward, and the band gap 
increases slightly from 0.34 to 0.38 eV upon the inclusion of the 
many-body effects. The slight increase in the band gap is attributed 
to the fact that the AM-covered silicene is a metal in fact, and the 
opened gap in the Dirac cone is 0.5 ~ 0.7 eV lower than Ef When Ef 
is shifted to the band gap region by the gate (off-state), the increasing 
magnitude of the band gap upon the inclusion of the many-body 
effects should be much larger. In the device simulation, the uncon- 
sidered many-body effects may cancel the artificial increment error 
of the transport gap by using the smaller SZ basis set^^. 

We also investigate the transport properties of the FET composing of 
a KSi2 or LiSi2 monolayer channel (Supplementary Fig. S6). When the 
SZ basis set is used, the transport gaps in the KSi2 and LiSi2 monolayer 
FET are calculated to be 0.7 and 0.6 eV, respectively, which are com- 
parable with that in the NaSi2 monolayer FET. The calculated on/off 
current ratio of the KSi2 monolayer FET is also up to 6 X 10^, suggestive 
of great potential of a KSi2 monolayer in a high-performance FET. 

In summary, our ab initio calculations demonstrate that AM sur- 
face adsorption is able to create a sizable band gap (> 0.4 eV) in 
silicene without degrading its extremely high mobility. Different 
from the band gap oscillating with the coverage in single layer gra- 
phene, the band gap opened in AM-covered silicene is dual-pro- 
tected by the breaking of sublattice and bond symmetry and 
generally increases with the coverage. The high on/off ratio of 10^ 
obtained in the simulated AM-covered silicene FET suggests poten- 
tial advantage of the AM-covered silicene for low-power consump- 
tion and high-speed logic devices. We anticipate experimental 
realization of the AM surface adsorption on silicene and outstanding 
performance in AM-covered silicene FETs. 

Methods 

First, we perform DFT calculations within the generalized gradient approximation 
(GGA) to the Perdew-Wang (PW91) exchange- correlation functional*". The geometry 
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optimizations and electronic structure calculations are performed with the all-electron 
double numerical atomic orbital plus polarization (DNP) basis sef\ as implemented in 
the DmoP package^^. A DFT-D semi-empirical correction is applied with the PW91 
functional to account for the dispersion interaction; the dipole correction is used to 
eliminate the artificial dipole moments interaction between the periodic images along 
the z-direction. The Monkhorst-Pack /c-point mesh^^ is sampled with a separation of 
about 0.01 in the BriUouin zone. The geometry optimization is carried out until the 
maximum force on each atom is less than 10"^ eV/A. In order to include the SOC and 
many-body effects of AM adsorption, we perform some further tests (For details, see 
Supplementary method). In a further step, the transport properties of the gated two- 
probe model is established by the DFT coupled with the nonequilibrium Green's 
function (NEGF) method, as implemented in the ATK 11.8 package^^'^l We employ 
the single-zeta (SZ) basis set during the device simulation, and a test with single-zeta 
plus polarization (SZP) basis set is also carried out for comparison. The Monkhorst- 
Pack /c-point meshes for the central region and electrodes are sampled with 1X50X1 
and 1X50X50 separately. The temperature is set to 300 K. Following the Landauer- 
Biittiker formalism^^, the current I is obtained by 

+ CO 

^(^g'^b) = ^ I {T{E,V„V^)\fiE-fi^)-f{E-^^)]}dE, (6) 

where T (E, Vg, Vb) is the transmission probability at a given gate voltage Vg and bias 
voltage Vb,/£) is the Fermi-Dirac distribution function, and fiJfiR is the electrochemical 
potentials of the left/right electrode (/^r — fii^ = eVb )• The electrostatic responses of 
and Vg are calculated by solving the Poisson equation self-consistently. The Neumann 
condition is used on the boundaries of the direction vertical to the silicene plane. On the 
surfaces connecting the electrodes and the central region, we employ Dirichlet boundary 
condition to ensure the charge neutrality in the source and the drain region. 
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